function r0 = set_init(initang)
% r0 = set_init(initang)
% This function sets the normalized coordintes to the corresponding angle

N=size(initang,2); % number of masses
rad=pi*initang'/180.;
srad=cumsum(rad);
r0=cumsum([cos(srad) sin(srad)])';

end